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1 Motivation 

Proper control of the numerical dissipation/filter to accurately resolve all rele- 
vant multiscales of complex flow problems while still maintaining nonlinear sta- 
bility and efficiency for long-time numerical integrations poses a great challenge 
to the design of numerical methods. The required type and amount of numerical 
dissipation/filter are not only physical problem dependent, but also vary from 
one flow region to another. This is particularly true for unsteady high-speed 

shock/shear/boundary-layer /turbulence/ acoustics interactions and/or combus- 
tion problems since the dynamics of the nonlinear effect of these flows are not 
well-understood. Even with extensive grid refinement, it is of paramount im- 
portance to have proper control on the type and amount of numerical dissipa- 
tion/filter in regions where it is needed. 


2 Objective 

The objective of this paper is to propose an integrated approach for the control 
of the numerical dissipation in high order schemes. Among other design criteria, 
the key idea consists of automatic detection of different flow features as distinct 
sensors to signal the appropriate type and amount of numerical dissipation/filter 
where needed and leave the rest of the region free of numerical dissipation con- 
tamination. These scheme-independent detectors are capable of distinguishing 
shocks/shears, flame sheets, turbulent fluctuations and spurious high-frequency 
oscillations. In addition, these sensors are readily available as an improvement 
over existing grid adaptation indicators. The detection algorithm is based on 
an artificial compression method (ACM) [2] (for shocks/shears), and redundant 
multiresolution wavelets (WAV) [1] (for all types of flow feature). The ACM 
and wavelet filter schemes using a second-order nonlinear filter with sixth-order 
spatial central scheme for both the inviscid and viscous flux derivatives and a 
fourth-order Runge-Kutta method are denoted by ACM66-RK4 and WAT 66- 
RK4. The scheme without any filter is denoted by CEN66-RK4 and the filter 
scheme using the dissipative portion of the fifth order WENO scheme (WEN05) 
as the nonlinear filter is denoted by WAVweno66-RK4. Computations blending a 
second-order nonlinear filter with an eighth-order linear filter using the wavelet 
as the sensor are denoted by WAV66-RK4-D8. For smooth flows, one of the de- 
sign criteria to improve nonlinear stability is the use of the entropy splitting of 
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the inviscid flux derivatives [3]. These schemes are indicated by appending the 
letters “ENT”. 

3 Numerical Examples 

Extensive grid convergence studies using WAV66-RK4 and ACM66-RK4 for a 
complex viscous shock/shear/boundary-layer interaction, and a viscous super- 
sonic combustion with four species (a planar Mach 2 shock in air interacting 
with a circular zone of hydrogen bubble) were conducted (see e.g., Fig. 1). 
More accurate solutions were obtained with WAV66-RK4 and ACM66-RK4 than 
with WEN05-RK4, which is nearly three times more expensive. Figure 2 shows 
an inviscid vortex convection problem at 150 periods on a 81 x 81 grid. The 
second-order MUSCL and the Tam & Web DRP method are extremely diffusive 
whereas WEN05 is more diffusive than the rest of the filter schemes using RK4. 
The CEN66-ENT-RK4-D8 linear filter scheme preserves the vortex convection 
beautifully. Figure 3 shows the highly accurate computation by blending of two 
different filters WAV66-RK4-D8 for a 1-D shock- turbulence problem. 

For a comprehensive study of this work, see [4]. 


WENO 1000x500 



(a) WEN05-RK4, 1000 x 500 


WAV 86 2000x1000 



Density at time 1 



(b) WAV66-RK4, 1000 x 500 


WAV68 4000x2000 



(c) WAV66-RK4, 2000 x 1000 (d) WAV66-RK4, 4000 x 2000 

Fig. 1 . Density contours of viscous shock/shear /boundary-layer interaction, Re 
1000 . 
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(a) Initial data 
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(d) CEN66-ENT-D8 
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(b) MUSCL 


ACM66-6NTk»0 01,150 P*t 81x81 po«nlx dl-0 4 



(e) ACM66-ENT 


VYAVWEt*O68-ENT|C-0 1t*»0.150 P*f BIxBl po*r»1» chO 4 



(h)WAVweno66-ENT 


Fig. 2 Vortex convection at 150 Periods, cfl— 0.4 (At 
MUSCL k DRP are extremely diffusive. 


DRP. d-0 02. 150 P*f 81x81 poxits ctoxO.l 
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(c) DRP 


WAV66-ENTi«0 01a»0, 150 P*f 81x81 poxitx.dWO 4 



(f) WAV66-ENT 


0.01; DRP 4 times smaller). 
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Solution, BSplines 



(a) 2nd-order UNO scheme on a 8000 (b) Density & wavelet detector (Al- 

grid (solid), WAV66-RK4-D8 on a 800 pha), WAV66-RK4-D8, 400 grid 

grid (circle) 

Fig. 3. Blending of two different filters by WAV66-RK4-D8 for a 1-D shock-turbulence 
problem. 
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